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The QCD ground state is highly non-trivial, and the strong condensates of scalar 
quark- ant iquark pairs (qq) and gluon fields (G c f lv G a,lv ) may play important roles in a wide 
range of low-energy hadronic phenomena 0, || . Therefore, it is quite interesting to 
study the density dependence of the condensates in nuclear matter. The vacuum values 
of the lowest-dimensional quark and gluon condensates are typically given by @ 

Q = (gg) ~-(225±25MeV) 3 , (1) 

G = (G%G afMU ) ^ (360 ± 20MeV) 4 . (2) 

Drukarev et al. ||, Cohen et al. || and Lutz et al. [BJ have shown that the leading 
dependence on the nuclear density, pb, of the quark condensate in nuclear matter, Q(ps), 
is given by the model-independent form: 

Q(Pb) , ctjv ,„v 

Qo fimi 

where <jn is the pion-nucleon sigma term (empirically <jn — 45 MeV |7j]), m % is the pion 
mass (138 MeV) and f„ ~ 93 MeV, the pion decay constant. Further, the strange quark 
content in the nucleon at finite density (and temperature) was studied in Ref. f| using 
the Nambu-Jona-Lasinio (NJL) model, supplemented by an instanton induced interac- 
tion involving the in-medium quark condensates. The gluon condensate at finite density, 
G(pb), has also been discussed in Ref. ||. 

Several years ago M. Ericson || suggested that a "distortion factor", coming from 
rescattering of soft pions in the nuclear medium, would tend reduce the amount of the 
chiral symmetry restoration (i.e., to oppose the reduction of the quark condensate). How- 
ever, Birse et al. [K| pointed out that there was an incompleteness in the treatment of 
the rescattering of soft pions using the simple linear a model, and showed that the full 
amplitude for soft-pion scattering from two nucleons leads to an enhancement of the chiral 
symmetry restoration at finite density. 

Recently the in-medium quark condensate has been calculated in several, more elab- 
orate ways [jO], |T2|. In particular, using the Dirac-Brueckner-Hartree-Fock (DBHF) ap- 
proach, Li and Ko and Brockmann and Weise have shown that higher-order contributions 



become increasingly important at high density, and tend to hinder the restoration of chiral 



symmetry 1 1 



We study the density dependence of the quark and gluon condensates in nuclear matter 



within the framework of the quark- meson coupling (QMC) model fl3| , [14], [T^, [16], [17]. The 
QMC model may be viewed as an extension of QHD |18| in which the nucleons still interact 
through the exchange of scalar (a) and vector (u and p) mesons. However, the mesons 
couple not to point-like nucleons but to confined quarks (in the nucleon bag). In studies 
of infinite nuclear matter it was found that the extra degrees of freedom provided by the 
internal structure of the nucleon give an acceptable value for the incompressibility once 
the coupling constants are chosen to reproduce the correct saturation energy and density 
for symmetric nuclear matter. This is a significant improvement on QHD at the same 
level of sophistication (see also Ref. [19[]). Furthermore, the model has been successfully 
applied to finite nuclei within the Born-Oppenheimer approximation [15 , |16| , [20| , ^TJ. It 
has been found that the QMC model can reproduce the properties of finite, closed-shell 
nuclei (from 12 C to 208 Pb) quite well. 

As shown in Refs. |To] , [RjJ, the basic result in the QMC model with mean-field ap- 
proximation (MFA) is that, in the scalar and vector meson fields, the nucleon behaves as 
if it had an effective mass M^. The latter can be calculated using a relativistic quark 
model for the nucleon (e.g., the MIT bag model) and depends on the nuclear density only 
through the a field. 



In an earlier version of the QMC model |15], |20| , we considered the effect of the nuclear 
medium on the structure of the nucleon alone and froze the quark degrees of freedom in the 
mesons. We call this version QMC-I. We have calculated the quark condensate in nuclear 
matter using this version [[14], | 



[23fl . (For a recent study, see also Ref. |[24| .) However, 
strictly speaking, those calculations (as well as the DBHF calculations fill) were not 
complete because the meson structure effects were not treated consistently. The mesons 
themselves are built of quarks and anti-quarks, and their structure may also change in 
matter [|I6"|, Ell. An additional, technical difference from Refs. IKI |22|, |23| is that for 
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reasons explained in Ref. ]T5[ - see especially Appendix - the cm. correction to the bag 



energy is now treated as being independent of the applied scalar field. 

To incorporate the effect of meson structure in the QMC model in MFA, we suppose 
that the vector mesons are again described by the MIT bag model with common scalar 
and vector mean-fields (like the nucleon in QMC-I). In this case the effective vector-meson 
mass in matter, m* (v = ui,p), will also depend on the a mean-field. The a meson itself 
is, however, not so readily represented by a simple quark model (like a bag), because it 
couples strongly to two pions and a direct treatment of chiral symmetry in medium is 



important 10, G5j. On the other hand, many approaches, including the NJL model 



the Walecka model |26| and Brown-Rho scaling |27| suggest that the u-meson mass in 
medium, m*, should be less than the free value, m a (= 550 MeV). It has therefore been 
parametrized as a quadratic function of the scalar field [f[6j| : 



a a {g*o-) + b a (g a cr) 2 , (4) 



with g a a in MeV. Here a is the mean-field value of the a field and g a is the a-nucleon 
coupling constant (in free space). Three parameter sets: (a a ; b a ) = (3.0, 5.0 and 7.5 xl0~ 4 
MeV" 1 ; 10, 5 and 10 xl0~ 7 MeV -2 ), called A, B and C, respectively, were determined in 
Ref. [IB | so as to reduce the mass of the a-meson by about 2%, 7% and 10% (respectively) 



at saturation density po (= 0.15 fm ). This model, involving the structure effects of both 



the nucleon and the mesons, was called QMC-II JT6 



Within QMC-II, the total energy per nucleon, E to t, can be written as |16] 



3« = —^Tz E dk^Mf + P + ^ + -4jp B + 7r4^pl (5) 
pB^Trf^J 2p B 2m* 2 8mfp B 

where g v {v = u>,p) is the f-nucleon coupling constant and kp t {% = proton or neutron) 
is the Fermi momentum. The density of protons (neutrons) p p (p n ) is defined by pi = 
kp./(3iT 2 ), and then Pb = P P + Pn and ps = p p — p n - Detailed values of the coupling 
constants and properties of nuclear matter in QMC-II can be found in Ref. (TB| . 

The density-dependent quark condensate, Q(pb), is formally derived by applying the 
Hellmann-Feynman theorem to the chiral-symmetry-breaking quark mass term of the 



total Hamiltonian. One finds the relation for the quark condensate in nuclear matter at 
the baryon density ps'- 



m q (Q(p B ) - Qo) 



d 



dm, 



-£(p 



B 



(6) 



where £(pb) = PBE to t and m q is the average, current quark mass of the u and d quarks. 
The resulting value for Q/Qo as a function of density is shown in Fig. [I] - dashed line. 
(We have chosen the quark mass to be m q = 5 MeV and the bag radius of the free nucleon 
Rn = 0.8 fm, but the result is quite insensitive to these choices.) 



Using Eq.(pf), the Gell-Mann-Oakes-Renner relation [28] and the explicit expression 



for the self-consistency condition of the a field in nuclear matter ||16| , Eq.fl5|) leads to the 
following explicit relation for the ratio of Q(pb) to Qq\ 



Q{Pb 

Qo 



1 



9a 



l 



6S N (0)mlfl 



P'r 



Pa 



9a \dm q 



(9aO-) 



III 



*2 



dm„ 



+ 



Po 



Am* 2 



(2/ P 



' dg 2 p 
dm n 



(7) 



where p r = ps/po, f P = P p /pb, Sjv(0) is the quark scalar charge of the free nucleon 
(= Jbagdf'^q'^q) and g% = g a / (3Sn(0)) (1^, |16|. The vector-meson mass is calculated 
using the bag model, while the cr-meson mass is given by Eq.(p. Because m q enters only 
in the combination m q — g^a, which is generally regarded as the chiral-symmetry-breaking 
term in nuclear medium, we were able to evaluate (^7-) in terms of the derivative of m* 
with respect to the applied scalar field a. 

In Eq.(^|) we have followed the usual convention of identifying 3m q 5 , 7v(0), which is the 
sigma commutator in the free MIT bag, as the experimental pion-nucleon sigma term, 
(Tat [|TJ|]. It is well known that the meson cloud of the nucleon (mainly the pions), as well 
as its strange quark content contribute significantly to [BIJ. However, because we are 
concerned primarily with the variation of Q in matter from its free value, Qo, it should 
be reasonable to replace <7at in Eq.(^) by its empirical value. (We note that the main 
variation of Q in medium is generated by the a mean-field.) 

Clearly, from Eq.(^), the leading dependence of the quark condensate on the density 
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is given by the scalar field: 

Q(Pb) 1 cr N fm c 



One can easily show that Eq.(Q) reduces to the model-independent result, Eq.(|3|), to 
leading order in the density, so that for small p r one has (for the set B) J2TJ : 

Q(p B )/Qo ^ 1 - 0.357p r . (9) 

This is shown as the dotted line in Fig. [l|. 

Equation (^) also involves deviations of the quark-meson coupling constants with 
respect to m q . In principle, if one could derive these coupling constants from QCD, 
their dependence on m q would be given. Within the present model there is no reason 
to believe that the couplings should vary with m q . This is especially so for the vector 
couplings since they involve conserved vector currents. On the other hand, we require 
that our model reproduces the correct saturation energy and density of nuclear matter 
whatever parameters are chosen for the free nucleon. As a consequence, the coupling 
constants depend on m q in a way that has nothing to do with chiral symmetry breaking. 
(For example, for set B, we find g q a = 4.891 - 0.005880m g + 1.200 x 10' 5 m 2 qJ gl = 
39.59 + 0.03828m 9 + 1.144 x 10 _3 m* and g* = 66.3 - 0.02m g , with m q in MeV.) 

In order to extract a physically meaningful result for Q/Qo we should therefore re- 
move the spurious contributions associated with ^ 1 (M = a,p,uS) in Eq.(0). In fact, 
the variation of g q with m q is extremely small so we need only correct the u and p con- 
tributions. The final, corrected result is shown as the solid line in Fig.|l| Even in the 
mean-field approximation, our calculations show that the higher-order contributions in 
the nuclear density become very important and that they weaken the chiral symmetry 



restoration at high density (c.f. Ref. [pT|). In QMC-II, the a field in nuclear matter 
is suppressed at high density (for example, g a a ~ 200 (300) MeV at p (3p )) because 
the quark scalar charge, Sn, decreases significantly as the density rises, as a result of the 
change in the quark structure of the bound nucleon [|TJ, |T(| |3D|. Since the reduction of the 
quark condensate is mainly controlled by the scalar field, it is much smaller than in the 
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Figure 1: Quark condensate at finite density using parameter set B. The dashed and 
dotted curves are respectively for the full calculation in symmetric nuclear matter (f p = 
0.5) and the linear approximation, Eq. ([)]). The solid curve is the result corrected by 
removing the spurious u) and p contributions. 



simple, linear approximation, Eq.([5|). From the difference between the solid and dashed 
curves we see that the correction for the dependence of the coupling constants on m q is 
significant and this should be born in mind in any phenomenological treatment. 

We should note here that, from extensive studies of chiral perturbation theory for 



nuclear matter, especially the recent work of Birse |L0], p5j , a reduction of the quark 
condensate from its vacuum value may not be enough to conclude that the chiral symmetry 
has been partially restored - especially if part of the change in (qq) arises from low- 
momentum pions. We note also that higher-order condensates may play an increasingly 
important role as the quark condensate tends to zero. 

Next let us consider the in-medium gluon condensate. Cohen et al. || also developed 
a model-independent prediction of the gluon condensate that is valid to first order in the 
nuclear density through an application of the trace anomaly and the Hellmann-Feynman 
theorem. 
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Following their approach, the ratio of the gluon condensate in nuclear matter, G(pb), 
to that in vacuum (G ) is given by 



G(p B )/G 



9G 



\£{Pb) - 2m q (Q(p B ) - Qq) - m s (Q s (p B ) - Q s o)) 



(10) 



where m s is the strange-quark mass and Q s (pb) (Qso) is the strange-quark condensate in 
nuclear matter (in vacuum). Up to first order in the density, the change of the strange- 
quark condensate may be written in terms of the strange quark content of the nucleon in 
free space, S: 

m s (Q s ( PB ) - Q s0 ) = Sp B + 0(p 2 B ). (11) 

The strange quark content is commonly specified by the dimensionless quantity, y, defined 
by 

ss Wh (12) 



(uu + dd) 



N 



which leads to S = (m s /2m q )<j]yy. Roughly speaking, y represents the probability to find 



s or s in the nucleon and is a measure of the OZI-rule violation. If m s /m q ~ 25 [28] 
and y ~ 0.45 ||, we get S ~ 250 MeV. We note, however, that y ~ 0.45 is an extreme 



value, and it has recently been suggested that it may be compatible with zero [[31]. In the 
analysis below we shall take care to examine the sensitivity to the full range of variation 
of y. 



At very low p B , £(pb) can be expanded as [|18 

3 



£(p B ) = M N p B 



1 + 



10M^ V 2 



+ 0(p 



B)i 



(13) 



where the second term in the bracket is the nonrelativistic Fermi-gas energy. Using the 
approximate form, g a a ~ 214 (MeV) xp r (for the parameter set B in m*) []2T|, we find 



2m g (Q(p B ) - Qo) = 214(MeV) x a N 1-^1 p r + 0(pf. 



(14) 



Choosing the central value of Go i n Eq-(@)j we then get the in-medium gluon condensate 
at low pb (for the set B): 



G(p B )/G = 1 - (0.03892p r + 0.001292p^ /3 ) + 0{p\ 



(15) 
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Figure 2: Gluon condensate at finite density using parameter set B. The solid and dotted 
curves are respectively for the full calculation in symmetric nuclear matter (f p = 0.5) and 
the approximation, Eq.flT5|). 

Our numerical results for the full calculation as well as the approximate calculation 
with Eq.(|T^) are shown in Fig. (for m q = 5 MeV and = 0.8 fm). The reduction of 
the gluon condensate at finite density is not large; for example, it is reduced by only 4% at 
Po, which is consistent with the results of Cohen et al. ||. The approximation of Eq.flTop 
works very well for a wide range of the nuclear density, which may imply that the effect 
of higher-order contributions (in powers of the density) is small for the gluon condensate. 
However, one should keep in mind that the in-medium gluon condensate evaluated here 
contains a large uncertainty, originating from the uncertainty in the value for the strange 
quark content of the nucleon in free space. We note that if we assume a vanishing strange 
quark content of the nucleon in free space (S = or y = 0) |3l| , the gluon condensate 
would be reduced by about 6% at po- 

Finally, we relate the quark condensate to the variation of the hadron mass in nuclear 
matter. In QMC-II, the hadron mass at low p# is simply given in terms of the scalar 
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field Pll 



M; ~ Mj - "j^cr), (16) 



where n is the number of non-strange quarks in the hadron j(= N, u,p,A, etc.). Since 
the quark condensate at low p_B is also determined by the scalar field (see Eq.(§)), we find 
a simple relation between the variations of the hadron mass and the quark condensate: 



r<fi\/g,\ 2 L Q(pb)\ , n L Q(p 



m *> ) "° I 1 " W ) " 2M(MeV) X "° I 1 ft ) ' (17) 



where 5Mf = Mj - M* (cf. Ref. [g7 



However, as shown in Ref. po[ , we know that the nucleon mass in matter cannot 
depend in any simple way on the quark condensate alone because the leading non-analytic 
contribution (LNAC) to the pion-nucleon sigma term - the term of order - should 
not appear in the nucleon-nucleon interaction To discuss this problem further, we 

have to include pions self-consistently in the QMC model, which is beyond the scope of 
the present work. 

In summary, we have calculated the quark and gluon condensates in nuclear matter. 
In the QMC-II model, the quark condensate at finite density is given in terms of the scalar 
field in the medium and the variation of the coupling constants with respect to the quark 
mass. We have shown that the reduction of the quark condensate at high density is much 
less than that suggested by the (model-independent) leading-order prediction, Eq.(|3]), even 
in the mean-field approximation. We also point out that the need to correct a naive use 
of the Hellmann-Feynman theorem to calculate Q{pb) for any purely phenomenological 
dependence of the quark- meson coupling constants on m q . 

In comparison with the quark condensate, the gluon condensate does not decrease 
much in nuclear matter. We have also provided a simple relationship between the change 
of the quark condensate and that of the hadron mass in nuclear matter. We should notice 
here that the validity of our model is limited to low and moderate density (probably less 
than ~ 3po), because the short-range correlations between quarks in overlapping hadron 
bags have been ignored. The effect of the pion cloud of the hadrons JlO|, [2R p3| should 
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also be considered explicitly in any truly quantitative study of the condensate properties 
in the medium. 
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